AUTP/04-07 



Stationary Metastability in an Exact Non-Mean Field Calculation for a Model 

without Long-Range Interactions 
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Abstract 

We introduce the concept of stationary metastable states (SMS's) in the presence of another more stable state. The stationary 
nature allows us to study SMS's by using a restricted partition function formalism as advocated by Penrose and Lebowitz and 
requires continuing the free energy. The formalism ensures that SMS free energy satisfies the requirement of thermodynamic 
stability everywhere including T = 0, but need not represent a pysically observable metastable state over the range where 
the entropy under continuation becomes negative. We consider a 1-dimensional m-component axis-spin model involving only 
nearest-neighbor interactions, which is solved exactly. The high-temperature expansion of the model representys a polymer 
problem in which m acts as the activity of a loop formation. We follow deGennes and trerat m as a real variable. A 
thermodynamic phase transition occurs in the model for m < 1. The analytic continuation of the high-temperature disordered 
phase free energy below the transition represents the free energy of the metastable state. The calculation shows that the notion 
of SMS is not necessaily a consequence of only mean-field analysis or requires long-range interactions. 



I. INTRODUCTION 



Supercooled and superheated states are ubiqutious in 
Nature, even though they cannot be rigorously derived 
from equilibrium statistical mechanics jl| . Their observa- 
tion is usually justified by appeal to the "van der Waals 
loop" in the celebrated van der Waals equation for the 
liquid-gas transition. The existence of the loop violates 
the fundamental property that the partition function 
(PF) be maximized or the free energy be convex. De- 
spite this, metastable states appear not only in many 
other mean-field theories such as the Bragg- Williams 
theory 0, but can easily be prepared in the labratory. 
There are usually two different mechnism operative in 
metastable states. There is usually a "fast" mechanism 
(time scale Tf) to create a metastable state in the system, 
followed by a "slow" mechanism (time scale r s ) for nucle- 
ation of the stable phase and the eventual decay of the 
metastable state. For the metastable state to exists for 
a while, we need to require r s > Tf. In approximate the- 
ories, the thermodynamic functions for the metastable 
states are taken as the extrapolation of the functions 
from the nearby equilibrium states. However, metastable 
states in real systems always have time-dependence asso- 
ciated with them. Thus, the metastable state represented 
by extrapolation can only represent the stationary limit 
(t s — ► oo)of experimentally observed metastable states; 
see also However, it has been suggested that the 
extrapolation is possible only because of the mean-field 
approximation, and would not be possible in real systems 
due to a singularity in the thermodynamic functions 0. 
The presence of the stable phase above some critical size 
in the metastable state is responsible for the decay of 
metastable states and for the essential singularity in the 
free energy The singularity is absent in mean-field 
theories or theories with long-range interactions. Ac- 
cordingly, the existence of SMS (no time-dependence) is 
commonly considered a mean-field consequence or due 
to long-range interactions so that one should not seen 



SMS's in real systems 0. On the other hand, what one 
observes in experiments are (time-dependent) metastable 
states. Since essential singularities are almost impossibe 
to detect experimentally, it is not surprising that the ex- 
trapolation is possible, at least from the experimentalist's 
point of view. 

In many cases, metastable states like supercooled liq- 
uids and glasses can remain stable for a long period of 
times 5, 6J. This should be contrasted with metastability 
at high temperatures in the liquid-gas transition that do 
not share this property. Thus, for supercooled liquids, 
t s >> Tf. This can be undestood by the high viscosity 
observed in supercooled liquids, which slows down the 
growth of the stable phase nuclei. There is another re- 
markable difference. Supercooled liquids such as viscous 
liquids usually do not (but very well could, as was seen re- 
cently 7]) exhibit spinodals, while supercooled vapor and 
superheated liquid invariably do. Rather, viscous liquids 
undergo a glass transition at low temperatures, about 
two-thirds of their melting temperature Tm, provided the 
liquid is cooled in a way that crystallization does not 
intervene. Here, the crystal phase (CR) represents the 
more stable phase, and care must be exercised to forbid 
its nuclei to form while cooling the viscous liquid. This 
makes the decay of the metastable state even less prob- 
able, and strengthens the inequality r s >> Tf. Thus, it 
is safe to treat viscous liquids as stationary metastable 
states (SMS's), which can then be described by equi- 
librium thermodynamics under the restriction that the 
crystal phase is not allowed. It is these SMS's that are 
of interest in this work. It is the hope that the study of 
SMS will throw some light on the properties of observed 
metastable states in the form of viscous fluids. In par- 
ticular, the extrapolated free energy below the melting 
temperature can be used to describe supercooled liquids. 

However, even if extrapolation is possible, one must 
still argue that the thermodynamic functions describe the 
stationary limit of experimentally observed metastable 
states. Under what condition(s) can one demonstrate 
this association to be valid? 
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FIG. 1: Schematic form of the generic entropy functions for 
various possible states. 



While there is no rigorous theory of such SMS's at 
present, there are some valuable approaches available in 
the literature. One such approach to describe SMS is 
to use the PL formalism of Penrose and Lebowitz (PL) 
1] using restricted ensemble method, which we modify 
and adapt for our case below. The modification is the 
following. The decay of the metastable states (to the 
stable state) in the PL formalism will be completely sup- 
pressed in order to make them stationary. Thus, nucle- 
ation of the stable phase will not be allowed in our study. 
This is consistent with Maxwell's idea [8| that to observe 
metastable states, we must ensure that the stable phase 
is not present. The properties of the SMS are what PL 
call the static or reversible properties 

An alternative scenario for extrapolation is by analyti- 
cally continuing the eigenvalues of the transfer matrix as 
presented in [!j , which attempts to accomplish the same 
as the restricted ensemble does but in a somewhat direct 
fashion. 



A. Schematic Entropy Functions under Continu- 
ation 

In the PL approach, only certain microstates out of 
all are allowed, the prescription of which is discussed in 
0. The restricted microstates are used to define a re- 
stricted partition function, which is then used to study 
metastable states. This is schematically shown in Fig. 1, 
where the curve OHAB represents the entropy function 
Soid(E) for the ordered crystal state, while DH AO K 
represents the entropy Sdi s (E) associated with the dis- 
ordered liquid state. The entropy as a function of E 
must be thought of as the entropy in the microcanon- 



ical ensemble [lOj, which must be at its maximum in 
the equilibrium state. Since a SMS is not an equilib- 
rium state in the unrestricted ensemble, its entropy at 
some E cannot exceed the entropy of the correspond- 
ing equilibrium state at the same E. It is clear, there- 
fore, that at lower energies, the ordered state must have 
higher entropy, while at higher energies the disordered 
state must have higher entropy. On the other hand, if 
a time-dependent metastable state is prepared under the 
constraint that the stable phase is not allowed, then the 
entropy function of such a metastable state will be rep- 
resented schematically by FG. The three free energies 
corresponding to the above entropy functions are shown 
in the inset. A consequence of the entropy maximiza- 
tion principle noted above is that the free energy F^i s (T) 
of SMS cannot be lower than the free energy F or( j(T) of 
CR at the same temperature T. This explains the form 
of the free energy in the inset. The slope of the tangent 
line HH gives the inverse melting temperature, while the 
slope of the tangent line 00 gives the inverse tempera- 
ture at which the free energy DO CK in the inset is equal 
to the free energy of the crystal phase at absolute zero 
(T = 0). The slope of the entropy at K is shown to be 
finite, as opposed to the infinite slope at O. This point 
will be discussed further below. 



The question that naturally arises is whether the above 
extrapolation is possible. It should also be noted that the 
extrapolation of the free energy does not guarantee that 
metastable states associated with this extension exist in 
the model. This will become clear in the following. Thus, 
the other important issue is to understand the condition 
under which the extrapolation will represent the station- 
ary limit of the metastable states that might be observed. 
To answer these questions, we borrow ideas from both 
approaches mentioned above and develop an approach, 
which is then tested by considering a 1-dimensional lat- 
tice model. This model has only nearest-neighbor inter- 
actions, and is solved exactly by the use of the trans- 
fer matrix. We find that the extrapolation can be car- 
ried out without any ambiguity to describe stationary 
metastable states (SMS) in this case. Thus, stationary 
metastability can exist even in non-mean-field theories 
and without long-range interactions, which is our main 
result. We further show that the extrapolation yields 
a thermodynamically stable SMS free energy, at least 
mathematically (see below for details) , all the way down 
to absolute zero. However, the continuation cannot repre- 
sent any metastable state at very low temperatures when 
the entropy becomes negative, and must be stopped. At 
this point, the continuation must be replaced by what 
is conventionally called an ideal glass; see below. This 
situation should be contrasted with the termination of 
a metastable state in a spinodal. The point where the 
entropy vanishes is not a spinodal. 
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B. Fundamental Postulate 



II. EQUILIBRIUM FORMULATION 



We assume the existence of SMS's, so that the partition 
function (PF) formalism can be applied. The need for the 
assumption is easy to understand. At present, our un- 
derstanding of whether equilibrium (lowest free energy) 
states can be demonstrated to exist mathematically even 
in simple models is too limited. We should recall that 
the existence of equilibrium states is taken for granted as 
a postulate in statistical mechanics and thermodynam- 
ics, where it is well known that it is extremely hard to 
prove their existence. We quote Huang "Statistical 
mechanics, however, does not describe how a system ap- 
proaches equilibrium, nor does it determine whether a 
system can ever be found to be in equilibrium. It merely 
states what the equilibrium situation is for a given sys- 
tem." Ruelle |l2j notes that equilibrium states are de- 
fined operationally by assuming that the state of an iso- 
lated system tends to an equilibrium state as time tends 
to +00. Whether a real system actually approaches this 
state cannot be answered. 

The problem becomes more complicated for SMS's like 
supercooled liquids in which, at least at low temper- 
atures, the relaxation becomes very sluggish and it is 
highly likely that the appropriate relaxation time T{ in- 
deed tends to +00. In other words, such an SMS may 
not even be observed in a finite amount of time, even 
though associated time-dependent metastable states can 
certainly be observed. Even in this case, the study of 
the long-time limit of metastable states still has a pre- 
dictive value, and can be carried out using the statistical 
mechanical formalism. 



C. Reality Condition 

For the microstates to exist in Nature, it is evident that 
W(E), the number of microstates of energy E, must sat- 
isfy the reality condition W(E) > 1 [so that the entropy 
S(T) > 0] even in the restricted ensemble. However, a 
state with negative entropy can emerge under extrapola- 
tions of the free energy. If it happens that the extrapola- 
tion results in a negative S(T) at low temperatures, this 
will indicate that the extrapolation no longer represents 
real microstates, and the system could not be found in 
those microstates in Nature. 

There are two independent aspects of thermodynam- 
ics and statistical mechanics. The first one is the re- 
quirement of stability according to which thermodynamic 
quantities like the heat capacity, the compressibility, etc. 
must never be negative. The other aspect, independent of 
the stability criteria, is the reality condition that ensures 
that such states occur in Nature [l^. The mathemati- 
cal extension of the free energy of the disordered phase, 
while always satisfying the stability criteria everywhere 
(T > 0), need not satisfy the reality condition, as our 
example will show below. 



A. Canonical Partition Function 

We consider a system composed of N particles confined 
in a given volume V and at a given temperature T. The 
canonical PF is given by 

Z N (T) ee Tr W N (E) exp(-0E), (1) 

where Tr is over all possible values of the e nerg y, Wn(E) 
is the number of microstates of energy E |14| , and /? = 
1/T,T being the system temperature in the units of the 
Boltzmann constant fee- We do not explicitly show the 
volume-dependence. We also introduce the adimensional 
free energy (without the conventional minus sign) 

n N (T) = lnZ N . (2) 

For microstates to exist in Nature, Wn(E) > 1; hence 
the corresponding entropy Sn(E) = In Wn(E) > 0. 
Whether this remain true for the analytic continuation 
remains to be seen. 



B. Thermodynamic Limit 

The thermodynamic limit is obtained by taking N — ► 
00, and V — > 00, keeping v = V/N fixed. The limit is 
taken by considering the sequence formed by 

u N (T) ee (1/N)Q N (T), 

for different values of N as N — > 00. The volume must 
be changed according to V = vN. For proper thermo- 
dynamics, the limit of the sequences must exist, which 
we assume and denote it by ui(T). The corresponding 
Helmholtz free energy is f(T) = -Tui(T). 

In the following, we will usually suppress the index TV 
on various quantities, unless necessary. 



C. Conditions for Equilibrium and Negative En- 
tropy 

We assume the existence of an equilibrium crystal, 
which has its energy E = Eq at T = 0. It also has 
the lowest free energy at low temperatures. Since Eq 
is an allowed energy, we must surely have W(E ) 7^ 0. 
Assuming TS(T) — ► as T ~ * 0, which is always true ac- 
cording to the Nernst's postulate, we recognize that Eq 
represents not only the Helmholtz free energy but also 
the energy of the perfect CR at T = 0. Since W(E) is 
non-negative, Z is a sum of positive terms. As a con- 
sequence, the following two principles of equilibrium are 
always satisfied. 
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1. Principles of Equilibrium 



III. STATIONARY METASTABLE STATES AND 
RESTRICTED ENSEMBLE 



• Maximization Principle The PF Z must be maxi- 
mized in the thermodynamic limit. The maximum 
value of Z(T) corresponds to picking out the maxi- 
mum term e s ~ pE in 0. This maximum term cor- 
responds to E = E. 

• Stability Principle The heat capacity, which is 
given by the fluctuations in the energy is non- 
negative. 

It should be stressed that the non-negativity of the 
heat capacity and the maximization principle only re- 
quire the positivity of W(E) (> 0); W(E) > 1 is not 
required. Thus, both principles remain valid even if the 
entropy becomes negative [13 ■ The above principles 
of equilibrium and reality are two independent aspects. 
This observation is going to be useful when we discuss 
the metastable states below. 



2. Principle of Reality 

Conventional statistical mechanics for a system in the 
thermodynamic limit describes equilibrium states in Na- 
ture, for which the above two principles of equilibrium, 
along with the principle of reality (W(E) > 1, S(E) > 0) 
must be satisfied All these conditions may not be 

met by metastable states. (Metastability does not occur 
in finite systems.) What we will see that it is the reality 
condition that can be violated by metastable states. 



D. Order Parameter 

The presence of a melting transition at Tm (the inverse 
of the slope of HH in Fig. 1) means that the disordered 
equilibrium liquid (EL) phase above Tm and the ordered 
CR below Tm correspond to different values of the order 
parameter p, which is traditionally defined in such a way 
that p = represents the disordered phase and p ^ 
the ordered phase CR. (Our example below will show 
explicitly how the microstates can be divided into the two 
disjoint classes.) We denote the free energy per particle 
above T M by Wdis(T) [ /dis(T) = -Tw d i S (T)], and below 
Tm by w or d(T) [ f ald (T) = -Tw ord (T)], from which we 
can calculate the entropies, and energies per particle 

s a (T) = -(dfc/dT), e a (T) = -(dw a /d/3), (3) 

a =dis, ord, respectively, corresponding to the two states. 
From s a (T) and e a (T), we can construct the functions 
s a (e) = s a [e a (T)], where e = E/N in the thermody- 
namic limit. The extensive entropy functions (s a multi- 
plied by N) are shown schematically in Fig. 1. 



A. PL Scheme 

We briefly review the restricted ensemble formalism 
developed by Penrose and Lebowitz pj , and the required 
modification to suit our purpose. Let ecR,M, and cel,m 
denote the energies of the coexisting phases CR and EL at 
the melting temperatute Tm; see points H and H in Fig.l. 
It is clear that Sdi s (e) and s or d(e) constructed above cer- 
tainly exist for e > eEL.M, and e < ecR,M, respectively. 
Over this range, we do not need to introduce the resticted 
ensembles. To obtain s a {e) beyond their respective range 
noted above, however, we need to introduce the restricted 
ensembles 1]. 

We begin by considering the case of finite but very 
large N. From s a (e), we can determine the number 
of microstates Wdi&(E) — exp[Nsdi s (E/N)] > 1 con- 
sistent with p = at high temperatures (or energies 
E > Tel.m = NeEL,M), and the number of microstates 
W OTd {E) '= exp[Ns rd(E/N)] > 1 consistent with p ^ 
at low temperatures (or energies E < Eqr,m = NecR,M)- 
(The equalities and inequalities are defined upto ther- 
modynamically insignificant terms.) Let us focus on 
Wdis(E) for E > Tel.m, which contains only those mi- 
crostates that are disordererd and correspond to p = 0. 
These microstates may contain a small number of clus- 
ters or nuclei of stable phase (CR), but their sizes are 
limited by the correlation length, which remains finite 
since we are dealing with a first-order transition. Let 
£dis (hi the units of some average inter-particle distance) 
denote the maximum value of the correlation length in 
the disordered phase. We now follow PL, and select all 
distinct microstates of energies E < Tel,m, in which 
there are no nuclei of the stable phase of sizes larger 
than £dis, and the number of smaller clusters is not too 
large, i.e. is thermodynamically insignificant to ensure 
that these configurations also correspond to p = 0; the 
check of this will be discussed below. We denote the num- 
ber of these microstates also by Wdi s {E). We can simi- 
larly extend W OI d(E) to E > Tcr,m- Thus, we can con- 
struct the two entropy functions S or d(E) = In Word (E), 
and Sdis(E) = \nWdi s {E) that overlap, and are shown 
schematically in Fig.l. 

B. Required Extension 

Let Eo denote the lowest energy in the system, which 
represents the energy of the ordered phase at T = 0. 
Thus, W otd {E) > 1. While W old (E) certainly exists for 
microstate energies starting from E = Eq, there is no 
guarantee that Wdis(E) also exists near E = E n . Most 
probably, Wdis(E) does not continue all the way down to 
E = E . If it did, the energy of the disordered phase at 
absolute zero would be E (we assume that TSdis — > 
as T — > 0), the same as that of CR. This would most 
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certainly imply that they would coexist at T = 0, each 
having the same volume; recall that we are considering 
a fixed volume ensemble. While there is no thermody- 
namic argument against it, it does not seem to be the 
case normally. Usually, the most stable state at T = is 
that of a crystal. Moreover, it is an experimental fact Q 
that all glasses have much higher energies or enthalpies 
compared to their crystalline forms at low temperatures. 
Thus, we assume that the lowest possible energy Ek , see 
Fig. 1, for the disordered state is larger than Eg. In other 
words, the microstate number W<n s (E) has the following 
property: 

W dis (E) > 1 for E>E K . (4a) 

If the slope in Fig. 1 at K is finite, then there is no 
sigularity in S,a s (E) at K, and we can extend it to lower 
energies. We assume this extension is possible and de- 
fine the extended entropy function for E > Eq. We de- 
note this extended entropy function by S dis (E), and in- 
troduce W£ is (E) = exp[Sl is (E)], The function S dis (£) is 
identical to Sdi s (E) over E > Ek- It exists over the entire 
range E > Eq, whereas Sdi s (E) exists only over the range 
E > Ek- We can similarly extend W or d(E) to E = Ej, 
where Ej is either equal to i?Max, the maximum allowed 
energy in the system, or the location of the singularity 
in S OI d(E) so that the latter cannot be extended beyond 
it. We denote this extension similarly by W* rd (E). In 
the following, we are mostly interested in the extension 

C. Restricted and Extended Restricted PF's 

Using W OI d{E), Wdis(E), and their extended version 
W* ld (E), W dis (E) we introduce the following restricted 
ensemble PF's Q: 

Z a (T) = TrW a (E) exp(-0E), (5a) 
Z* a (T) = TrW* a {E) exp(-/?£), (5b) 

a =dis, ord, and the corresponding free energies 

n a (T) = ]nZ a {T), n* a (T) = lnZ*JT). (6) 

The free energy per particle f2 Q (T) /N is expected to pos- 
sess a thermodynamic limit as N — * oo, which we have 
already introduced earlier as oj a (T). The corresponding 
limiting free energy per particle O* (T)/N will be denoted 

by<(T). 

• Remark The following remark is important to un- 
derstand the relationship between the starred and 
unstarred PF's. Let us consider the disordered 
PF's. For temperatures so that the average energies 
E dis (T) and E dls (T) are greater than Ek, both par- 
tition functions are determined by the microstates 
of energies above Ek, where the starred and un- 
starred W dis s are identical. Hence, for T > 2k, the 



two PF's Z dis (T) and Z dis (T) are the same so that 
their free energies are the same. They differ only 
below T K ; while Z dis (T) exists there, Z dis (T) does 
not. Similarly, for temperatures so that E ord (T) 
and E old (T) less than Ej, Z* md (T) and Z ord (T) are 
the same. Thus, 

u dis {T) = u* dis (T), T>T K , (7a) 
w OTd (T)=u4 d (T), T < Tj. (7b) 

Here, Tj is the temperature where E OTd (T)— Ej. 
The free energies iv dis (T) is defined for all temper- 
atures T > 0. 

As long as W a (E) > 0, and W*(E) > 0, the restricted 
PF's are sum of positive terms. Therefore, the corre- 
sponding free energies satisfy the two equilibrium condi- 
tions noted above. Consequently, even the restricted and 
extended restricted PF's will never give rise to unstable 
states. 

It is clear that the global maximization of the PF re- 
quires that 

uj(T) = cjdis(r) = u dis {T), T > Tm, 
w(T) = w ord (T) = w* rd (T), T <T M . 

The switchover from cj d i S (T) to w ord (T) at Tm makes 
oj(T) singular, as expected, due to the transition. 

We consider the case when there is only one phase 
transition, the first-order melting transition, in the sys- 
tem. The following point is to be noted as discussed by 
Penrose and Lebowitz Q]. The resticted PF's defined in 
(|5aH and (|5b|) require that we add an extra energy term 
in the energy of the system, which takes the value if 
the microstate belongs to the set a, and +oo, if it does 
not. This meets the PL criterion for "static" metastable 
states. The other two criteria that PL require relate to 
the decay of metastable states, and does not have to be 
imposed here anymore. Thus, the problem of two incom- 
patible requirements discussed by Penrose and Lebowitz 
1] no longer is an issue. 

A prescription to describe metastability using the PF 
formalism can now be formulated. 



D. Metastability Prescription 

We abandon the above global maximization princi- 
ple, and use w d ; s (T) to give the free energy of the 
metastable disordered phase (supercooled liquid) below 
Tm and u> ord (T) to give the metastable (superheated crys- 
tal) state free energy above Tm- Similarly, s d j s (T), e d ; s (T) 
and s OT d(T), e d ; s (T) give the entropy and energy per par- 
ticle for the supercooled liquid and superheated crystal, 
respectively. 

There are two possibilities for the extrapolation of the 
free energy. As said above, unstable states are not pos- 
sible in the restricted ensemble. Thus, either the free 
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energy terminates in a spinodal at a non-zero but finite 
temperature, or it extrapolates to T — through u dis (T) 
for the supercooled liquid (T — > oo through w* rd (T) for 
the superheated crystal) . In this work, we are only inter- 
ested in the supercooled liquid. 

It is easy to calculate the order parameter p for T < Tm 
for the "disordered phase" by using Z^ is (T) to check if we 
have properly identified the set of disordered microstates 
above. Since all microstates in W d i s (E) contain only nu- 
clei of the stable CR phase of finite sizes, the argument of 
Fisher Q about the origin of an essential singularity no 
longer works, which requires nuclei of all sizes, including 
infinitely large sizes. Thus, it is clear that lo^IT) can 
be used to describe the saught extrapolation of the free 
energy below the melting temperature. The phase rep- 
resented by u) dis {T) below Tm will still correspond to a 
disordered state (p = ). This is our required descrip- 
tion of SMS in the form of SCL by the PF Z dis (T) below 
Tm- 

From the above discussion, it appears highly likely that 
the singularity in uj(T) does not necessarily imply a sin- 
gularity in either of its two pieces oj d i s (T) and uj ord (T). 
Both of them can exist on either side of Tm- From the 
above argument, we conclude that the extrapolation used 
to define Z* a (T) is not a consequence of any approxima- 
tion (mean- field or otherwise). Our example below is 
intended to give a concrete demonstration. 

The form of the entropy functions S a (E) shown in Fig. 
1 is also supported by all known observations 5, ,6J , exact 
calculations 0, flfij l , from the arguments given above 
and the calculation to be presented below. We note that 

S* old (E) < S dis (E), E>E M , (8a) 
S* dis (E) < S ord (E), E < Em, (8b) 

where Em is the energy at A where S md (E) = S d i s (E); 
see Fig. 1 . The SMS corresponding to the stationary 
SCL is defined by the branch H ACK and its extention 
to Eq, which is not shown. Similarly, superheated CR is 
defined by the branch HAB and its extension to higher 
energies. We note that, as shown, the entropy S d [ s of the 
metastable branch goes to zero at Tk>0 corresponding 
to the finite slope at K. This behavior will be supported 
by the exact calculation in the next section. 

IV. EXACT 1-D CALCULATION 

The calculation presented here follow the transfer ma- 
trix eigenvalue approach of Newman and Schulman 
We now consider a one-dimensional axis spin model, 
which contains m-component spins Si located at site 
i of the one-dimensional lattice of N sites, with peri- 
odic boundary condition (Sat + i =Si). Each spin can 
point along or against the axes (labeled 1 < k < m) 
of an rn-dimcnsional spin space and is of length : 
S = (0, 0, .., ±y/rn, 0, ..0). The spins interact via a ferro- 
magnetic nearest-neighbor interaction energy (— J), with 



FIG. 2: The bond and the entropy densities. The bond den- 
sity is a monotonic function of T, so that the stability is not 
violated. The entropy becomes negative at low temperatures, 
where metastable state must be replaced by an ideal glass. 



K = J/T > 0. The energy of the interaction is given by 
E = — J Si ■ Sj+i. 

i=l...N 

The PF is given by 

/ 1 \ N ( 1 \ N N 

Z N{ K,m)={-) EcxpR^-J Tr T , 

(9) 

where the first sum is over the (2m) N spin states of the 
N spins and T = exp(iifS • S') is the transfer matrix 
between two neighboring spins. The transfer matrix has 
the eigenvalues 

Adi S = u + 2(m - 1), A ord = v, A = u — 2, (10) 

that are 1-fold, m-fold, and (m— l)-fold, respectively [l8j . 
Here we have introduced the following 

x = exp(i^m), u = x + l/x,v = x — 1/x. 

We follow de Gennes 0,H(j and provide an alternative 
and very useful interpretation of the above spin model in 
terms of a polymer system, in which each polymer has 
multiple bonds and loops. The valence at each site in 
a polymer must be even. (The presence of a magnetic 
field will allow odd valencies, which we do not consider 
here.) The high-temperature expansion of the PF, which 
is given by 

Z N (K,m) = Y J K B m L , (11) 

describes such a polymer system, with K > 0, and m de- 
noting the activity of a bond and the activity for a loop, 
respectively, and B and L denoting the number of bonds 
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and the number of loops, respectively |20j. The empty 
sites represent solvent particles. The number of polymers 
and the number of bonds and loops in each polymer are 
not fixed and vary according to thermodynamics. In ad- 
dition, there is no interaction between polymers, and be- 
tween polymers and solvent particles, so that the polymer 
system in (| 1 1|> is an athermal solution. The temperature 
T of the spin system does not represent the temperature 
in the polymer problem, as is well known [T^. |20). We 
will see below that small x corresponds to high tempera- 
tures where the disordered phase is present, and large x 
corresponds to low temperatures where the ordered and 
possible SMS phases are present. Thus, decreasing T 
amounts to going towards the region where the ordered 
and metastable disordered phases are present. Let ui de- 
note the limiting value as N — > oo of 



lo n = {l/N)\nZ N {K,m) +ln(2m), 



(12) 



where we have added an uninteresting constant to get rid 
of the prefactor in This is done because the number 
of microstates appears within the summation in the spin 
model PF in 0. Thus, the inclusion of the prefactor will 
make the microstate entropy negative. The prefactor is, 
however, required for the polymer mapping. 

The importance of the polymer mapping is that we can 
take m > to be a real number, even though non-integer 
m makes no sense for a physical spin. Thus, for non- 
integer values of m, only the polymer system represents 
a physical system. For m = 1, the axis model reduces 
to the Ising model, while for m — > 0, it reduces to the a 
model of linear chains with no loops [T^. l20j . The eigen- 
value Adis is dominant at high temperatures for all m > 
and describes the disordered phase. Its eigenvector is 

(Xdisl = 

where (2k\ ( or (2k + 1|) denotes the single-spin state in 
which the spin points along the positive (or negative) fc-th 
spin-axis. It has the correct symmetry to give zero mag- 
netization (p = 0). For m > 1, Adis remains the dominant 
eigenvalue at all temperatures T > 0. For < m < 1, 
the situation changes and A or d becomes dominant at low 
temperatures T < T c , or [x > x c — 1/(1 — m)] where T c is 
determined by the critical value x c = exp(Jm/T c ); there 
is a phase transition at T c . The corresponding eigenvec- 
tors are given by the combinations 



(fe+i) 

Aord 



= [(2k\ - (2fc + l|]/V2, 



(k = 0, 2, .., m — 1) which are orthogonal to (xdisl, as can 
be easily checked. These eigenvectors have the symmetry 
to ensure p ^ 0. The remaining eigenvalue A is (m — 1)- 
fold degenerate with eigenvectors 



(fe+i) _ 



= [(2k\ + (2k + 1| - ((2k + 2| - (2k + 3|)]/V4, 



(k = 0, 2, .., m — 2.) For m > 0, this eigenvalue is never 
dominant. For m — > 0, it becomes degenerate with Adis- 



Since the degeneracy plays no role in the thermodynamic 
limit, there is no need to consider this eigenvalue sepa- 
rately for m > 0. 

We now consider the limit N — > oo. The adimensional 
free energy per site, which represents the osmotic pres- 
sure 0, l2l| , of the high-temperature equilibrium phase is 
^dis(T) = ln(Adis)- It can be continued all the way down 
to T = 0, even though the equilibrium osmotic pressure 
has a singularity at x c . Similarly, u OT d{T) = ln(A or d) re- 
lated to the low-temperature equilibrium phase can be 
continued all the way up to T — > oo. To calculate the 
entropy density, we proceed as follows. The bond and 
loop densities are given by 



6b = doj/d\n K, 0l = duj/dlnm, 



(13) 



which are needed to calculate the entropy per site of the 
polymer system 

s^ F ' = uj — (f)B In K — </>l In m; 

the superscript is to indicate that it is the polymer system 
entropy, and is different from the spin system entropy 
S (S) _ QTuj/dT. If we define lj without the last term in 
(tT2"l) . then </>l and must be replaced by (</>l — 1) and 
(s( p ) — In 2), respectively. This will not affect any of the 
conclusions below. 

In the following, we will be only interested in the poly- 
mer entropy. The proper stability requirements for the 
polymer system are 



(d^/dlnK) > 0, (d(j> h /dlnm) > 0, 



(14) 



as can easily be seen from (JHJ), and must be satisfied even 
for SMS. They replace the positivity of the heat capacity 
of the spin system, which no longer represents a physical 
spin system for < m < 1. It is easy to see from the 
definition of that (ds^J /dT) m need not be positive, 
even if the conditions in (|14f> are satisfied. 

Let us compute uj as K — > oo (T — > 0) for the two 
eigenvalues Adis and A or d- From l|13|) . it is easy to see 
that <t>B ~ * rnK for both states as T — > 0. Thus, using 



S (P) 



lnK 



, In m, we have 



w d is(T)/w ord (T) — > 1 as T 



0. 



(15) 



This means that if the eigenvalue Adis is taken to repre- 
sent the metastable phase above x Cl its osmotic pressure 
must become equal to that of the equilibrium phase (de- 
scribed by the eigenvalue A or d) at absolute zero. This is 
in conformiity with Theorem 3 in We take u),a B (T) 
to represent the SMS osmotic pressure below T c . We have 
also checked that Tsf^, -> 0, as T -> 0. 

We will only discuss the disordered polymer phase be- 
low for < m < 1. It is easily checked that the above sta- 
bility conditions in (|14fl are always satisfied for Adis! see, 
for example, the behavior of <^b in Fig. 2, where we have 
taken m = 0.7, and ./ = 1. Since the high-temperature 
disordered phase represents a physical system, it cannot 
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give rise to a negative entropy s dis above T c ; however, 
its metastable extension violates the reality principle as 
shown in Fig. 2, where its entropy s^ is becomes nega- 
tive below Tk — 0.266, which is lower than the transition 
temperature T c . 

We now make an important observation. As m de- 
creases (below 1), both Tk and T c (Tk < T c ) move down 
towards zero simultaneously. As m — > 0, the equilib- 
rium ordered phase corresponding to A or d disappear com- 
pletely, and the disordered phase corresponding to X dls 
becomes the equilibrium phase. There is no transition 
to any other state. Thus, there is no metastability any- 
more. Consequently, there is no ideal glass transition 
since there is no other state more ordered than this state 
any more, as argued above. Thus, our exact calculation 
confirms our earlier conclusion that the existence of an 
ordered state is crucial for the existence of the entropy 
crisis. The existence of an ordered state sets the zero of 
the temperature scale by its minimum energy E . This 
scale then sets the temperature Tk of the lowest SMS 
energy TJk > Eq to be positive. 

We also observe that there is no singularity in Adi s or 
^dis(T) at T c , even though there is a phase transition 
there. Similarly, there is no singularity in A or d or w or d(T) 
at T c . Thus, the thermodynamic singularity in the equi- 
librium free energy does not necessarily create a singular- 
ity in a;dis(T) or w or d(T) at T c , as was discussed earlier. 
The existence of a singularity or spinodal at some other 
temperature is a different matter. 

V. DISCUSSION & CONCLUSIONS 
A. SMS & Exact Calculations 

The transition between SMS and the ideal glass is not 
brought about by any thermodynamic singularity at Tk; 
rather, it is imposed by the reality requirement. The ideal 
glass state does not explicitly emerge as a new phase in 
the calculation since it is a disordered phase defined by 
the order parameter p = 0. In this sense, the transition to 
the ideal glass is a very special kind of transition, which 
does not seem to belong to the class of phase transitions 
in which various phases emerge in the calculation. 

The exact calculation, which is not mean-field calcula- 
tion in principle, in the previous section demonstrates the 
existence of SMS. Thus, it demonstrates that our hypoth- 
esis of SMS existence is not vacuous. It also shows that 
the free energy can be extrapolated below the melting 
temperature by the use of the restricted PF, and that 
there is no essential singularity, a signature of a first- 
order transition Q. The free energy remains stable all 
the way down to absolute zero. However, the mere exis- 
tence of the stable extrapolated free energy all the way 
down to T = does not mean that it represents the free 
energy of a realizable metastable state. This becomes 
evident when we consider the entropy of the extrapo- 
lated free energy. This entropy drops rapidly, and goes 



through zero at Tk, and becomes negative as the temper- 
ature is reduced. A genuine entropy crisis appears in the 
SMS below T K . At T K , /| iB (T K ) = E K , which is higher 
than the CR free energy Eg at T = 0. Below Tk, the 
extrapolated free energy cannot represent any real met- 
satable state and must be replaced by another free energy 
branch, which is constant: /ig(T) = Ek for T < Tk- It is 
shown by the dotted horizontal straight line at K in the 
inset. This branch represents the free energy of the ideal 
glass (IG), which is the phase below Tk- We need to in- 
voke an ideal glass transition at this temperature in the 
model. The energy of the ideal glass is Ek-This means 
that the ideal glass has a higher energy than the crystal 
at absolute zero, in conformity with the experiments. 

It is interesting to note that Tk — * 0, as m — ► 0, so 
that the ideal glass transition disappears. This is not 
surprising, as T c — > 0. Thus, the there is no ordered 
state anymore. 



B. No Entropy Crisis in the Equilibrium State 

The entropy crisis occurs only in the metastable state, 
and not in the equilibrium state, even though we have 
not shown this explicitly here. The entropy of the latter 
vanishes at Eq with an infinite slope, as shown in Fig. 1. 
Thus, the lowest energy Eq determines the lowest allowed 
temperature T = in the system, which is consistent 
with the Nernst-Planck postulate. However, it is possi- 
ble that the equilibrium free energy becomes horizontal, 
so that the entropy vanishes, over a non-zero temperature 
range (0, Tc) where the system is frozen. Such a situation 
happens, for example, in the KDP model and gives rise 
to a singularity at Tq, This should be contrasted with 
the existence of the ideal glass transition in the super- 
cooled state, where its appearance is not accompanied 
by any singularity in the SMS free energy. Replacing 
the unphysical SMS free energy below Tk by a frozen 
state is done by hand; it does not emerge as part of the 
calculation. Indeed, as our calculation has shown, the 
ideal glass transition disappears as m — > 0. In this limit, 
T c — * 0. Thus, the "ordered" state corresponding to A or d 
disappears, and the disordered phase remains the equilib- 
rium state all the way down to T = 0. Thus, it is safe to 
conclude that equilibrium state in any system will never 
show an entropy crisis (at a positive temperature). If 
any exact calculation for the free energy or the entropy 
predicts an entropy crisis at a positive temperature, this 
will necessarily imply that there must exist another state, 
the equilibrium state, for which no entropy crisis should 
exist. 

This observation has been crucial in a recent investiga- 
tion of a dimer model j2^| in which the disordered phase 
underwent a first-order transition to an equilibrium or- 
dered phase. The ordered phase then gave rise to an 
entropy crisis at a lower temperature, which forced us to 
look for another equilibrium state, which was eventually 
discovered above the temperature where the entropy cri- 
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sis was found, so that the crisis occurred in a metastable 
state (this time emerging form an intermediate ordered 
state). Similar situation occured in more complex sys- 
tems containing particles of different shapes and sizes 

In summary, we have shown that stationary metastable 



states can appear in exact calculations also. They do not 
only occur in mean-field calculations. 

We would like to thank Andrea Corsi and Fedor Seme- 
rianov for various useful discussions, and help with the 
first figure (Andrea Corsi). 
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